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Description of ^-distributions at RHIC energies in terms of the 
Ornstein-Uhlenbeck process and relativistic diffusion equations 
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1 Department of Physics, Shinshu University, Matsumoto 390-8621, Japan 
2 Toba National College of Maritime Technology, Toba 517-8501, Japan 
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Assuming the Ornstein-Uhlenbeck process and relativistic diffusion equations with three 
sources at y — 0, we analyze ^-distributions of d+Au collision as well as Au+Au collisions 
at ^/s]7n = 130 GeV and 200 GeV. 

1 Introduction: ^-distributions in Au+Au collisions at RHIC with various centrality- 
cuts (C.C.) at RHIC energies have been reported in Refs. 1) and 2). In Ref. 3) we 
have proposed a stochastic model with the sources at ± y m ax and y = 0, to describe 
the ^-distributions. That is an extended model from a model of two sources at 
±y,max in Refs. 4) and 5) In this letter, we propose a model with three sources 
at y = 0. A more realistically physical picture is shown in Figs. 1 and 2, where our 
basic assumption is shown therein. 

Very recently PHOBOS Collaboration has reported ^-distribution on dN/drj (d 
+ Au) at ^/I^a7 = 200 GeV. 6 ) Thus we are interested in analyses of the data, 
to examine whether or not the model proposed in Figs. 1 and 2 works well for 
descriptions of dN/ drj (d+ Au) as well as dN / dr\ (Au+Au) . 

2 Theoretical formula: To improve the physical picture of Ref. 3), we consider a 
formula based on the stochastic approaches. To describe physical picture of Figs. 1 
and 2(b), we have found that the following relativistic diffusion equations are avail- 
able from studies of Ref. 6), 



dP(y, t) 
dt 



d_ 

dy 



(y±y r , 



+ 



la 2 d 2 



2 7 dy 2 



P(y, t) 



Using the initial condition 

P(y, 0) = 5(y) , (2) 
we obtain the following solutions for the projectile and target fragmentations 



■. exp 



(y - ymax (i 



- 7t))2 



2V 2 (t) 



(3a) 



*' In Ref. 5), we has pointed out that single Gaussian distribution with an initial condition S(y) 
cannot explain the whole ^-distributions in Au+Au collisions. This fact suggests that contributions 
from the projectile and target fragmentations are necessary. 
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Fig. 1. Description of two-step processes by relativistic diffusion equations and the Ornstein- 
Uhlenbeck process. 




Fig. 2. Description of two-step processes for d+Au collisions. 



Umax j t) 



: exp 



2vm 



(3b) 



where V? p or T) = (a 2 /2-/)(l - e" 27 *). For the 3rd source due to gluons in the central 
region, we assume the Ornstein-Uhlenbeck process (Eq. (1) without the term y ma x', 
see Ref. 3)). The following solution can be derived from Eq. (2) as 



P {c) (y, t) 



exp 



2V \t) 



(4) 
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Combining Eqs. (3) and (4), we obtain the following expression 
tr, = J(r? ' m/Pt)X { c T+ C c P + c P(T){y > Vmax > f) 

+ J P - P (p) (y, -y max ,t)\ 

CT + Cp + Co J 

+Jo\ J\ PW(y,t)\, (5) 
[ C T + C P + c J 

where cp, cp and Co are the weight factors of three contributions for the target 
fragmentation, the projectile one and the central region, respectively. A relation 
between y and rj, and the Jacobians are given as 

I a/1 + (m/pt) 2 + sinh 2 rj + sinh 77 

y=2 ln ^— ' (6) 

1 + (m/pt) 2 + sinh 2 77 — sinh r/ 



a 1 \ cosh 7? 

J(r?, m/pt) = —. (7) 
y 1 + (m/p t ) 2 + sinh 2 77 

and J = J(r], m /p t ). 

3 redistributions in Au+Au collisions at smaller centrality-cuts (C.C.): 

We apply Eq. (5) to the data on 77-distributions at the C.C.6-15% and 15-25% at 
./sj^ = 130 GeV and 200 GeV by the PHOBOS Collaboration. Our results, using 
the CERN MINUIT program, are shown in Fig. 3 and Table I. We have to determine 
the weight factors, cp, cp and Co- For this aim, we adopt the following criterions 
m o/p< > 5( m o/Pt), rn/pt > S{m/p t ) and mo/pt < Sm/pt, because pions are mainly 
produced in the central region. They reveal that the ratio is cp : Co : cp = 1 : 6-7 : 1 
at the minimum x 2 's, where the evolution time p = (1 — e -27 *) stops. 



Table I. Values of estimated parameters of Eq. (5) with cp/cr = 1 for Au + Au at ^/snn = 130 
GeV and 200 GeV. 

130 GeV 200 GeV 



c. c. 


6-15% 


15-25% 


6-15% 


15-25% 


V 


0.52±0.21 


0.62±0.13 


0.78±0.12 


0.66±0.16 




5.17±0.95 


5.98±1.03 


5.64±1.42 


6.06±1.67 


vS 


1.31±0.26 


1.34±0.30 


1.43±0.25 


1.59±0.28 


CO 


2465±39 


1808±32 


3078±54 


2120±40 


cp 


c /6 


co/7 


co/7 


co/ 6 


mo/pt 


0.44±0.26 


0.37±0.29 


0.47±0.33 


0.42±0.33 


m/pt 


3.60±2.30 


2.34±1.06 


1.06±0.88 


2.06±1.14 


X 2 /d.o.f. 


0.64/48 


0.82/48 


0.50/48 


0.54/48 



4 redistributions in d+Au collision: Taking into account the mass ratio, : 
rriAu = 1 : 100, we can fix the ratio as cp/cp = 100. The weight factor co is treated as 
a free parameter in the present analysis. In eq. (5), two Jacobian factor, J(rj, rriT/pt) 
and J(rf, mp/p t ) are used. The evolution parameter p = 1 — e -27 * is also stopped at 
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minimum \ 2 - Our result is given in Fig. 4 and Table II. We know that the weight 
factors are described by cp : Co : ct = 0.53 : 40 : 53. This ratio means contribution 
from the gluon is the same magnitude of the target fragmentation. 



Table II. Parameter values of Eq. (5) with c T /c P = 100 for d+Au at ^/s^7 = 200 GeV. 6) 



p 


Vi 


v2 


v£ 


C() 


0.76 


4.27±0.64 


0.165±0.159 


7.56±1.38 


39.8±5.4 


cp 


mo/pt 


mp/pt 


mr/pt 


x7d-o.f. 


0.529±0.077 


0.556±0.166 


1.63±0.69 


1.88±0.41 


5.1/45 



5 Concluding remarks: 

1) To improve the previous physical picture in Ref. 3), we assume the relativis- 
tic diffusion equations for the target and projectile fragmentations in addition 
to the Ornstein-Uhlenbeck stochastic process available for that in the central 
region. 

2) The ratio, cp : cq : cp, is determined by means of the following criterions; 
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Fig. 4. Analyses of data on d+Au collision at ^/sjvjv = 200 GeV by Eq. (5). 

mo/pt > S(mo/pt), m/pt > S(m/pt) and mo/pt < Sm/pt- They can be investi- 
gated by analyses of 77-distributions with restrictions 

Vci ~| v\ in a future. 

3) As seen in Figs. 3 and 4, our formula seems to be available for the description 
of 77-distributions in d+Au as well as Au+Au collisions. 

4) From analysis of d+Au collision, cp : cq : ct = 0.5 : 40 : 53, we have know 
that the contribution of gluon in Au is the same ordered magnitude as that of 
quarks. 

5) Since our formula is expressed in the analytical form, it can be applied to other 
reactions, for example, Si+Al and S+Au collisions. These analyses will be 
shown elsewhere. 8 ^ 
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